New mechanism to cross the phantom divide 
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Recently, type la supernovae data appear to support a dark energy whose equation of state w 
crosses — 1, which is a much more amazing problem than the acceleration of the universe. We show 
that it is possible for the equation of state to cross the phantom divide by a scalar field in the gravity 
with an additional inverse power-law term of Ricci scalar in the Lagrangian. The necessary and 
sufficient condition for a universe in which the dark energy can cross the phantom divide is obtained. 
Some analytical solutions with m < -1 or tu > -1 arc obtained. A minimal coupled scalar with 
different potentials, including quadratic, cubic, quantic, exponential and logarithmic potentials are 
investigated via numerical methods, respectively. All these potentials lead to the crossing behavior. 
We show that it is a robust result which is hardly dependent on the concrete form of the potential 
of the scalar. 

PACS numbers: 95.36.+X 98.80.Qc 04.20.Cv 



I. INTRODUCTION 



Over the course of the past decade, the evidence for the acceleration of the the universe [l| is one of the most 
striking cosmological discoveries. A number of approaches, namely the existence of dark energy, have been adopted to 
. try and explain the remarkable observation. An incomplete list includes: quintessence, phantom, k-essence, modified 
Oh' gravity, and so on The quintessence field models correspond to an equation of state w > — 1. The region where 
the equation of state is less than — 1 is typically refer to as a presence due to something called phantom dark energy. 
Furthermore, in the wake of the more accurate data, a surprising result emerges: the resent study of the type la 
supernovas data show that the time varying dark energy provides a better fit than the cosmological constant, in 
addition, especially, the equation of state parameter w, defined as the ratio of pressure to energy density, may cross 
the phantom divide w = — 1. This crossing to the phantom region is neither possible for an ordinary minimally 
£T} ■ coupled scalar field nor for a phantom field. It clearly shows that this transition is difficult to be realized by only a 
^ ' single-field models of dark energy. In this article, we try to explore new ways and means to obtain the consistency for 
, the theoretical outcome with the present observation. 

The contribution of the matter content of the universe is represented by the energy momentum tensor on the right 
hand side of Einstein equations, whereas the left hand side is represented by pure geometry. Thus, there are two ways 
to give rise to current acceleration: (i) supplementing the energy momentum tensor by an exotic matter with negative 
pressure such as a cosmological constant or a scalar field; (ii) modifying the geometry itself. Since the dark energy 
models mentioned previously have not gotten fully satisfactory explanations of cosmic acceleration, it is worthwhile 
considering the possibility that acceleration of universe is not due to some kind of stuff, but rather is caused by 
modified gravitation. Lately, there is a very interesting operation on the general relativity: a modification to the 
Einstcin-Hilbert action by involving new terms of inverse powers of the curvature scalar (in the following text we call 
it inverse-R gravity), of the form gR~ n (with n > 0) [3(. Terms of this form would become important in the late 
universe and can bring about self-accelerating solutions, supplying a purely geometric candidate for dark energy. But, 
as shown in [3|, the effect of extra term yJ—gR~ n equals to a scalar field with a smooth positive potential. Therefore, 
the EOS of dark energy still fails to cross the phantom divide yet. Given the challenge of this problem, we investigate 
a novel model to pursue the challenge. In general f(R) gravity (even without scalar field), the phantom divide crossing 
behavior of dark energy is shown in 0, [B| , and i? 2 -term can yield a transient crossing [j, Q . 

In the next section we shall make a brief review of inverse-R gravity and its relation to the scalar theory in Einstein 
gravity. In the section III we investigate the mechanism for a single scalar to cross the phantom divide by an analytical 
method. In section IV, we present the numerical results of the crossing behavior of a scalar with different potentials 
in frame of inverse-R gravity in detail. And in Section V, we will present the main conclusion and some discussions. 



II. INVERSE-R GRAVITY 



There are varies modifications of Einstcin-Hilbert action. Basically, we classify them into two categories: one is 
ultraviolet modification which will be important in the high curvature region (early universe) while the other is infrared 
modification which plays significant role in the low curvature region (late universe). The ultraviolet modification has 
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been investigated widely in the inflation models, for example the starobinsky inflation [6|. Recently, the infrared 
modification arouses much interests because of the discovery of cosmic acceleration. The inverse-R gravity is one of 
the leading model in the model with infrared modifications. It is noted that the inverse-R gravity is equivalent to a 
class of scalar tensor theory with ui — 0, which is not compatible with solar system observations 0- However, the 
physical meaning of this approach is not clear enough as the first sight. For example, it is shown that a specific f(R) 
gravity implying cosmic acceleration does not physically equal to scalar tensor theory, though it do equal scalar tensor 
theory mathematically @. Thus, it is still sensible to study modified gravity and scalar tensor theory, respectively. 

Moreover, a scalar curvature squared term to the action can save the inverse-R gravity and help it to pass the 
observations in solar system Q . We study the behavior of a scalar field in the inverse-R gravity with R 2 correction, 
whose action S reads, 

S=^Jd* x V- 9 (f(R)) + I &xy/=g {Cm + C phi ). (1) 

Here m p i is Planck mass, g denotes the determinant of the metric, R marks the Ricci scalar, a, f3 are constants with 
mass dimension, Cm is the matter Lagrangian, 



f(R) = R- — +P 2 R\ (2) 



and Cfj, labels the Lagrangian of a scalar field with potential V ((/>), 



c^ = --d^d^-v(4>). (3) 



In this set-up, (3 is at the scale of inflation and a is at the scale of the present universe, hence j3 ^> a. Note that 
R 2 term aids the original inverse-R gravity to pass the tests in solar system and to ensure the universe undergoes a 
matter-dominated era, though we should ensure its effects in the late universe are tiny. The action ([lj yields field 
equation, 



where TM,,T^ V denote the energy- momentums of matter and scalar field, respectively. In the present article, we 
suppose that matter is just dust. By assuming an FRW universe with scale factor a, 

ds 2 = -dt 2 + a 2 dT, 2 , (5) 

where dT, 2 denotes the metric of a 3-dimensional maximally symmetric space, we obtain the friedmann equation 
corresponding to (UJ), 

2 k _ F(H, H, H) _ p M + 2 /2 + V 

a 2 3m 2 pl ~ 3m 2 pl ' W 

where H = a/a is the Hubble parameter, a dot denotes derivative with respect to the cosmic time t, and 

2 4 9 

F = ^ ( 2HH + 15H 2 H + 2H 2 + 6iJ 4 ] - 18(6H 2 H - H 2 + 2HH)-£-. (7) 

12(H + 2H 2 f V J 1 ' /3 2 w 

Another equation to close this system is the acceleration equation, which can be replaced by the equation of motion 
for <j), 

4> + 3H<p+^- = 0. (8) 

dep 

In this construction the scalar field together with the extra geometric term plays the role of dark energy. First we 
present a concise note on the definition of dark energy. In the inverse-R gravity theory, there is a surplus F-term 
in the modified Friedmann equation ©. However, almost all observed properties of dark energy are obtained in 
frame of general relativity. To explain the the observed evolving EOS of the effective dark energy, we introduce the 
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concept "equivalent dark energy" or "virtual dark energy" in the modified gravity models [10( . We derive the density 
of virtual dark energy caused by the non-minimal coupled scalar by comparing the modified Friedmann equation in 
the inverse-R gravity with the standard Friedmann equation in general relativity. The generic Friedmann equation in 
the 4-dimensional general relativity can be written as 

where the first term of RHS in the above equation represents the dust matter and the second term stands for the 
dark energy. Comparing © with ©, one obtains the density of virtual dark energy in STG, 

Pde =P<f>+F, (10) 

where 




We see that the F-term makes a key difference from the standard model. In the next section, we shall study the 
dynamics of a universe in inverse-R gravity, especially the dynamics of the virtual dark energy (dark energy for short). 



III. EXACT SOLUTIONS IN INVERSE-R COSMOLOGY 



To find exact solutions is an important but difficult topic in such a high non- linear system © and ([5]). In this 
section we study some exact solutions to explore the mechanism to cross the phantom divide for the dark energy in 
inverse-R cosmology. In this section we only consider an ideal case, that is, ft — 0, to see some quanlative behavior of 
the late universe in inverse-R gravity. And we leave the more realistic case with i? 2 -term in the numerical analysis. 



A. Power-law solution 



In a spatially flat universe, the Friedmann equation © reduces to 



2 k F(H, H, H) 2 /2 + V 

H + — o 2 = o 2 ) \ 12 ) 

a 2 3m z pl 3m' pl 

during the dark energy dominated stage. Square-law solution arises when the potential is chosen to take the expo- 
nential form, 

y = m 4 e 4l ; (13) 
where m e is a constant. With this potential, (fT2|) and (jHJ have a particular solution as follows, 

a = a Q t 2 , (14) 



<j) = 2m pl \n(^— ). (15) 
V 10 rripi ) 

It is easy to calculate that w = —2/3 in this case. One sees that this universe is "uniformally accelerating", since 
a =constant. An interesting point in this uniformally accelerating universe is that the extra term F vanishes just in 
time, and hence it is the solution with which we are familiar in standard general relativity by coincidence. 



B. De Sitter solution 



The dynamical system of the scalar filed with canonical and non-canonical Lagrangian has been widely studied 
among which the global structure of the phase space has been investigated and various critical points and their 
physical significance have been identified and manifested. We have shown that if the potential of the model with 
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kinetic energy (quintessence) has a non- vanishing minimum [12J, or the potential of a model with negative energy 
(phantom) has non- vanishing maximum (l3j , the dynamical system of the model admits late-time attractor solutions 
corresponding to — — 1 and 0^ = 1, where il denotes the fraction of a stuff of the universe. 
For a class of model that admits late-time de Sitter attractor, we have a solution, 

a = aoe [^/6™? i +(V;, 2 /r^ ! +a 4 /12) 1/2 /2] 1/2 (t-to) ; ( lg ) 



(17) 



From (fT6l) and (0, we derive 



a ml, 



Vo^ + ivi/gm^ + ayn)^ 



(18) 



Thus = 4^- = and the crossing phenomenon cannot happen near the late-time de Sitter attractor. 



C. Big rip solution 

Present observation data indicate the possibility of dark energy with w < — 1, dubbed as phantom [l4|. In the 
ordinary case, the realization of w < — 1 could not be achieved. Unfortunately, it suffers from the well known quantum 
stability problem. Though there are a few discussions to evade the stability problem for phantom model, we are far 
from "solving it" A good w < — 1 model should avoid the negative kinetic energy as much as possible. 

Under the potential, 

where I{u) is the inverse function of the function u = m p i((3\nv H — \— ), and 

/3 = 7 + V73, (20) 



we have a big rip solution, 



200(25 + V73) ' 
a = a (t~ t hr y 2 . 



(21) 



(22) 



4> = m p i 



3 7 (i - t br ) 4 



(23) 



where tb r is the rip epoch. For this solution, we have 



(24) 



D. The crossing mechanism 

Since there are both w > — 1 and w < — 1 solutions, the crossing phenomenon is bounded to happen in our model. 
Before studying the crossing behavior in this model, we say a bit more on the relation between density evolution 
and the equation of state. In fact, we can know that dark energy behaves as quintessence, vacuum or phantom only 
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from its rate of change of the density. We reach this point from the following arguments. Since the dust matter obeys 
the continuity equation and the Bianchi identity keeps valid, dark energy itself satisfies the continuity equation 



dpde 
dt 



+ 3H(p de +p ef f) = 0, 



(25) 



where p e f / denotes the effective pressure of the dark energy. And then we can express the equation of state for the 
dark energy as 



Peff 
Pde 



= -1 - 



dpd 



3pde da 



(26) 



From the above equation we find that the behavior of Wde is determined by the term — (cosmological 

constant) bounds phantom and quintessence. More intuitively, if pd e increases with the expansion of the universe, 
the dark energy behave as phantom; if p de decreases with the expansion of the universe, the dark energy behave as 
quintessence; if pde decreases and then increases, or increases and then decreases, we are certain that EOS of dark 
energy crosses phantom divide. A more important reason why we use the density to describe property of dark energy 
is that the density is more closely related to observables, hence is more tightly constrained for the same number of 
redshift bins used [161 ] . With data accumulation, observations which favor dynamical dark energy become more and 
accurate. Furthermore, a significant possibility appears recently: the EOS of dark energy may cross —1 (phantom 
divide) [13], which is a serious challenge for theoretical physics. The theoretical explore of the cros sing phenomenon 
was proposed in 18 1 . This interesting topic is under intensively studying very recently , see also [ldj ] and for some 
earlier references therein. 

can be rewritten as 



w = — 1 — 



1 



3Hp de 



■n 



(27) 



The behaviors of F and p^, increase or decrease with expansion of the universe, depend on the rate of expansion. A 
simple calculation presents, 



F 



Hm 2 pl a 4 



6(2H 2 

For a power-law universe, 



H) 



H(2H 2 + H)+ 3H(2H 2 - 7HH - H) + 3H(H 2 - 16H 2 H - 4ff 4 ) 



(28) 



do 



(29) 



where to is the present age of the universe, £ is a constant. £ > denotes an expanding universe; £ > 1 implies an 
accelerating one; £ > 2 implies a > 0, hence F > and F > 0, which presents a possibility to cross the phantom 
divide for the dark energy. Substituting (|2"9")l into (j2"B)l . we derive, 



F 



2m -i) 2 



In the £ = case, we have F max = 1.75 x 1CT 6 . If £ >> 1, we have F « 1. It is consistent that F 

de Sitter expansion, because de Sitter expansion is faster than any power-law accelerating expansion. 
From © and (H]), we find that 



(30) 



in the 



H 



-3H<f> 2 + F 



(31) 



H = 



-3H(f> 2 



6## + F — &m 2 pl H 2 



6m 2 pl H 



(32) 



The universal argument is that the crossing phenomenon happens at t = t c iff H t= t c = and H t= t c 7^ 0. In the 
quintessence case in standard model, ie, F = 0, H t=tr — implies H t=tc 7^ 0. Therefore, the crossing behavior never 
happens. In inverse-R gravity, Ht=t c = and H t=ta ^ iff, 



F = 3H<t> 2 , 



(33) 
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and 

F 3H<j) 2 . (34) 

We now have proved two theorems as follows, 

Theorem 1 The crossing phenomenon happens iff F = 3Hcb 2 , F ^ QH<fxf), <j> ^ and cj> ^ 0. 

Theorem 2 Iff F < 3H(f> 2 , we have the model w < — f ; iff F > 3H<f> 2 , we have the model with w > — 1. 

IV. NUMERICAL EXAMPLES 

In the above section, we discuss the evolution of the universe in the scalar dominated stage via a analytical way. 
For a realistic universe, the pressure-less dust is a necessary component. However, it difficult to study the evolution 
of the universe when we introduce a dust in inverse-R gravity with R 2 corrections. In this section, we investigate the 
dynamics of scalars with five different potentials, including quadratic, cubic, quantic, inverse power-law, exponential 
and logarithmic potentials in a universe with dust component via numerical methods, respectively. In the following 
text, we always assume a spatially flat universe. 

A. quadratic potential 

Quadratic potential is the most widely-investigated potential in field theory, which represents a mass term. Any 
potential around its minimum (if it has a minimum) can be treated as quadratic potential. In this subsection we 
study the case of quadratic potential, 

V = m 2 <j> 2 1 (35) 

where to is a parameter with dimension of mass. 

For convenience, we introduce two dimensionless variables 

x^H/H ,y = cj)/M p , (36) 

with which we rewrite Friedman equation and the equation of motion (JSJ) for ip into the following form, 

x 2 = IW- 3 * + \x 2 y 12 + \k 2 y 2 + h{Ax2x ' 2 X*f X " + 15 2 f 3 X ' + 6xi) + 18c(6xV + x 2 x' 2 + 2s V), (37) 
6 3 3b(xx' + 2x A Y 

= xx'y' + x 2 y" + 3x 2 y' + 2k 2 y, (38) 
where a prime denotes the derivative with respect to the so-called e- folds s = In a, b = c = is a dimensionless 

2 

constant, the coefficient k 2 = -ni and 



^ = -wFt- (39) 



Here Hq is the current Hubble parameter, pmo represents the current density of dust. Generally, c is a tiny number 
since the inflation scale f3 is much higher than Ho. Hence, the effects of i? 2 -term in ^ can be omitted in late universe, 
although it becomes important in the early universe and structure formation era. This also implies that the exact 
solutions discussed in the last section make physical sense in the late universe. 

Though we get the sufficient condition for the F-term to help the dark energy to cross the phantom divide, it is 
difficult to apply this condition before we obtain the exact solution. However in a universe with multiple components, 
one hardly obtain an exact solution. Under this situation, we present a numerical result about the dark energy density 
in Fig.l. For the sake of convenience, we introduce the dimensionless density as below, 



3iiQTOp Z 6 3 36{xx' + 2x A y 
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FIG. 1: /3, w and q as functions of s. This figure is for V = m 2 (f> 2 , in which we set Q m o = 0.3, = 1, 6 = 20, c = 10 116 , and 
initial values for arguments: a;o = 1, yo = 0,x' o — —0.3,y' = 0. 

The most important parameter from the aspect of observation is the deceleration parameter q, which carries the 
total effects of the cosmic fluids, 
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-1 



x 

X 



(41) 



We also displays q in Fig. ((T|). 

From Fig. ([IJ, obviously, the equation of state of dark energy crosses -1. Simultaneously, the deceleration parameter 
is consistent with observation. As is known to all, the equation of state of a single scalar in standard general relativity 
never crosses the phantom divide, the _F-term, induced from the term in the Lagrangian, plays an essential role 
in this crossing. 

B. cubic potential 

Next, we study the potential V = m^cf) 3 . This potential has been investigated in classical theory, especially in the 
dynamics of an oscillator. As for V — TO3<^ 3 , the Friedman equation §§§ and the equation of motion ([5]) become 



1 



x 2 = Q m0 + -x 2 y a + ^-k 3 y 
6 3 



1 3 b{4x 2 x' 2 + 2x i x" + 15arV + 6x 4 ) 



36(xx' + 2x 2 ) 3 
= xx'y + x 2 y" + 3x 2 y' + 3fc 3 y 2 , 



18c(6arV + x 2 x' 2 + 2x 3 x"), (42) 

(43) 



where &3 = rripimz/H 2 . 

The results in this situation is presented in Fig. 2. One can see easily the crossing behavior as mentioned above. 



C. quantic potential 



In this subsection we explore the dynamics of the universe for a scalar with V ~ ccj) (c is a dimensionless constant) 
in inverse-R gravity. The quantic potential is extremely important in modern physics. Due to Higgs mechanism, every 
massive particle get mass through a quantic potential. Then, for V = c<jr , the two principle equations become 



1 



1 



x 2 = f2 m0 + 7^V 2 + xfc42/ 4 
6 3 



b{4x 2 x' 2 + 2x 3 x" + 15x 3 x' + 6x 4 ) 
36(xx> + 2x 2 ) 3 



+ 18c(6xV + x 2 x' 2 + 2x 6 x") 



(44) 
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in which — crn^J Hq. 

We give the results in Fig. 3. 
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FIG. 4: This figure is for V = m\e m v l , f2 m o = 0.3, and initial values for coefficients: fc e = X, b = 10, xq = l,yo = 0,x' o — 
-0.3,y' o =0,c= 10" 118 . 

D. exponential potential 

In this subsection, we study the exponential potential. In the standard inflation models, the exponential potential is 
an important example which can be solved exactly in the standard model. In addition, we know that such exponential 

potentials of scalar fields occur naturally in some fundamental theories such as string/M theories. We present the 

t 

results about V = m^e™? 1 in following content. 

Concerning V = m^e m p' , the equations © and flSJ) to be 

x 2 = Q m0 + \x 2 y' 2 + \k e ey + K^ 2 x' 2 +2x^' + 15x 3 x' + 6x^ + ^ + ^ + ^ 
6 3 So(xx' + 2x*y 



= xx'y' + x 2 y" + 3x 2 y' + k e e v , 



(47) 



where k e — m^H m pl . 

The results are given by Fig. 4. 



E. logarithmic potential 

The logarithmic potential is not widely studied in particle physics. To fit parts of elementary particle mass spectra 
by involving logarithmic potentials comes into physicists' view only in recent decades [20]. In this paper, we study it 
phenomenologically. With regard to V = mf In -2- , the corresponding equations ([6]) and (HJ will be 



x 2 = tt m0 + \x 2 y 2 + -ki \ny 
6 3 



b(4x 2 x' 2 + 2x i x" + 15arV + 6x 4 ) 
36(xx> + 2x 2 ) 3 



+ 18c(6x d x' + x 2 x u + 2x A x"), (48) 



= xx'y' 



x 2 y" 



3xV + -- 



(49) 



where ki — mfH 2 m p 2 ■ 
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We present the results in Fig. 5. 

From figs. 1-5, one sees that all the w of the virtual dark energy with different potentials are successfully to make the 
equation of state of dark energy cross — 1. And the 3 curves in all the figures have a similar shape, which implies that 
the crossing behavior is determined by the extra geometric term. Consequently, we have to say that it is geometric 
property of the inverse-R gravity itself, which is independent of the individual potentials. 



V. CONCLUSION 



In conclusion, the results of this paper demonstrate convincingly that it is possible to realize the crossing w = — 1 
for the equation of state of the dark energy by a single scalar. We find the necessary and sufficient condition for a 
universe in which the dark energy cross the phantom divide in inverse-R gravity. And then wc investigated different 
potentials to minimally coupled scalar field <fr , including quadratic, cubic, quantic, exponential, logarithmic potentials 
in inverse-R gravity with R 2 correction. And the results state clearly that different potentials lead to the crossing 
behavior, respectively. Therefore, we conclude that it is a robust property of inverse-R gravity with R 2 correction, 
not controlled by a special potential. 
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